Quantum and semiclassical study of magnetic anti-dots 
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We study the energy level structure of two-dimensional charged particles in inhomogeneous mag- 
netic fields. In particular, for magnetic anti-dots the magnetic field is zero inside the dot and 
constant outside. Such a device can be fabricated with present-day technology. We present detailed 
semiclassical studies of such magnetic anti-dot systems and provide a comparison with exact quan- 
tum calculations. In the semiclassical approach we apply the Berry- Tabor formula for the density 
of states and the Borh-Sommerfeld quantization rules. In both cases we found good agreement 
with the exact spectrum in the weak magnetic field limit. The energy spectrum for a given missing 
flux quantum is classified in six possible classes of orbits and summarized in a so-called phase dia- 
gram. We also investigate the current flow patterns of different quantum states and show the clear 
correspondence with classical trajectories. 

PACS numbers: 73.21.-b, 03.65.Sq, 85.75.-d 



I. INTRODUCTION 

In the past decade, the study of systems of two- 
dimensional electron gas (2DEG) in semiconductors^ 
has been extended by the application of spatially inho- 
mogeneous magnetic fields. The inhomogeneity of the 
magnetic field can be realized experimentally either by 
varying the topography of the electron gaa^, by us- 
ing ferromagnetic material a'^1^1^1^1^1^ , depositing a super- 
conductor on top of the 2DECUS*ii. Numerous the- 
oretical works also show an increasing interest in the 
study of electr on motion in inh omogeneous magnetic field 
(see, e.g., Refs Jl2ll3ll4ll.5llfil lr7ill8.ia.2a.2L.22.23.24J25L 
m 27, 28 29 30il^ 

In the experimental works mentioned above, for GaAs 
heterostructures, on the one hand, the electron dynamics 
is confined to two-dimensions. On the other hand, the 
coherence length and the mean free path of the electron 
can be much larger than the size of the system, while the 
Fermi wavelength is comparable to the size of the 2DEG. 
Moreover, the electron system can be described to a good 
approximation as a free electron gas with an effective 
massi. Therefore, the quantum mechanical treatment of 
these systems is of some physical interest. 

In this paper, as an example, we consider the energy 
levels of a two-dimensional non-interacting electron gas 
in a magnetic field that is zero inside a circular region 
and constant outside. This system (shown in Fig.^ will 
be called a magnetic anti-dot; it was first studied by Soli- 
many and Krameri2^. Solving the Schrodinger equation it 
was shown that there are bound states. Introducing an 
effective angular momentum, the Schrodinger equation 
of the particle in symmetric gauge can be mapped to the 
Landau model. This effective angular momentum is a 
sum of the angular momentum in a uniform magnetic 
field and the flux (in units of the flux quantum) miss- 
ing from the uniform field. Recently, Sim et ali2^ have 
renewed the study of this system and pointed out the 



crucial role of the magnetic edge states in the magneto- 
conductance. The classical counterparts of these states 
correspond to trajectories of the charged particle that 
consist of straight segments inside the non-magnetic re- 
gion and arcs outside. 




FIG. 1: The two-dimensional electron gas in an inhomoge- 
neous magnetic field. The magnetic field is zero inside the 
circle of radius R and constant outside the circle. 

Although it is difficult to measure directly the den- 
sity of states of a quantum system, it affects many ob- 
servable quantities such as the magnetoconductance, the 
magnetization or the susceptibility. In the interpretation 
of the experimental results the semiclassical approxima- 
tion proved to be a useful tool. Several semiclassical ap- 
proache9Si2^iSi2Si2L22i2£ are known in the literature and 
an excellent overview of the subject can be found in the 
textbook by Brack and Bhaduri^. Different semiclassi- 
cal theories for magnetic systems have successfully been 
applied, for example, in works^ii^Siiiii. For integrable 
systems Berry and Tabor— have shown that the oscillat- 
ing part of the density of states can always be expressed 
in terms of classical periodic orbits. This formula is com- 
monly called the Berry- Tabor trace formula. 

One of our aims in this paper is to apply, for the 
first time, the Berry- Tabor trace formula for a magnetic 
anti-dot. To illustrate the power of the method, we 
also calculate the exact eigenvalues of the single parti- 



cle Schrodinger equation and find a very good agreement 
between the two results. We should mention here that 
the statement by Sim et ali^S on the relation between the 
quantum states and the corresponding periodic orbits is 
somewhat misleading. Their condition for a given peri- 
odic orbit is not necessarily satisfied at the value of the 
corresponding exact energy level as they claimed How- 
ever, including more and more periodic orbits with the 
proper weights in the trace formula the sum converges to 
the correct quantum density of states. In practice only 
a few of the shortest orbits are enough to get a rough 
estimate of the positions of the exact energy levels. 

The power of the semiclassical approach can also be 
demonstrated by applying the Bohr-Sommerfeld approx- 
imation. We shall show that the energy levels obtained 
from the Bohr-Sommerfeld quantization rules also agree 
very well with the numerically exact levels even for the 
lowest eigenstates. Note that in this semiclassical treat- 
ment the quantization should be applied to the classi- 
cal motion on a two-dimensional torus parametrized by 
the action variables and their canonically conjugate angle 
variables (for details see, eg, Ref. 40). 

The classical orbits can be classified by their cyclotron 
radius g and their guiding center c (distance of the center 
of the orbits from the origin). In the quantum mechani- 
cal treatment one can calculate the average of operators 
defining the cyclotron radius and the guiding center. For 
circular magnetic billiards, Len^^ has derived approxi- 
mate expressions for these averages for a given quantum 
state. Following Ref. l45l one may derive the correspond- 
ing relations for magnetic anti-dot systems. Thus, these 
relations are the basis for classifying the different quan- 
tum states in terms of classical orbits for our system. We 
will show that the quantum states can be described by six 
different types of classical orbits. In addition, this clas- 
sification enables us to draw a so-called 'phase diagram' 
which shows a clear one-to-one correspondence between 
classical orbits and quantum states. 

To complete our semiclassical study, we finally present 
results for the probability current density calculated from 
quantum calculations. We shall argue that the current 
fiow patterns can be understood qualitatively from the 
corresponding classical trajectories. Recently, Halperin^ 
has been shown that the total current (the integral of the 
current density along the radial direction) can be related 
to the dispersion of the energy levels (their angular mo- 
mentum dependence). As it will be shown this general 
relation works in our magnetic anti-dot system, too. 

Regarding the numerical calculations, we should men- 
tion that the semiclassical approach presented in this 
paper proves to be a very effective method. Moreover, 
it provides a better understanding of the nature of the 
quantum system. Our semiclassical method applied to 
magnetic anti-dot systems may be an important tool to 
understand the role of the magnetic edge states in the 
density of states or the magnetization (both are experi- 
mentally accessible physical qauntities). We believe that 
our semiclassical analysis can be extended to other types 



of inhomogeneous magnetic fields such as studied, eg, in 
Refs. .(LT- S^Q. as well as non-circular dot systems. 

The rest of the text is organized as follows. In Sec. 
^the exact quantization condition (secular equation) is 
derived from the matching conditions of the wave func- 
tions at the boundary of the magnetic and non-magnetic 
regions. In Sec. IIIII the semiclassical approximation is 
presented including the description of the classical mo- 
tion of the particle in Subsec. IIII Al the characterization 
of the possible periodic orbits in Subsec. 1111 Bl some nu- 
merical results in Subsec. IIII CI and the phase diagram 
in Subsec. IIIIDI The current flow patterns of the sys- 
tem are discussed in Sec. lIVI Finally, the conclusions are 
given in Sec.lVl 



II. QUANTUM CALCULATION 

In this section, we present the quantum mechanical 
treatment of the magnetic anti-dot. The magnetic field 
with a constant B outside a circle of radius R is as- 
sumed to be perpendicular to the plane of the 2DEG. 
The Hamiltonian of the electron of mass M and charge 
e is given by 



H 



p 

2A/' 



2M 



if r<i?, 
, if r > R, 



(1) 



where p is the canonically conjugate momentum, and the 
vector potential in polar coordinates (r, ip) and symmet- 
ric gauge is given by"^^ 



A = A^{r,ip)e^, where 
r"^ -R^ 



(2) 



and e^ is the unit vector in the (p direction. Here Q{x) 
is the Heaviside step function. 

The energy levels of the system are the eigenvalues E 
of the Schrodinger equation: 



H-^{r,Lp) = E^{r,ip). 



(3) 



Rotational symmetry of the system implies a sepa- 
ration ansatz for the wave function as a product of 
radial and angular parts. We choose for the angular 
part the appropriate angular momentum eigenfunctions 
^imip -^j^jj quantum number m (here m is an integer). 
Thus the wave function for a given m is separated as 
'i'{r,ip) — fm{''')e""''^ , where the radial wave functions 
fm{r) satisfy a one-dimensional Schrodinger equation in 
the normal region: 



Ki{T)fm{T) =ef,n{T), 

in which the radial Hamiltonian takes the form 
0)2 1 d 



hm{T) 



dr^ T dr 



+ V^{t). 



(4a) 



(4b) 



Here we introduce the dimensionless variable r 



r/U 

where I = \Jh/\eB\ is the magnetic length, Wc = \eB\/M 
is the cyclotron frequency, e = 2E/(hujc) is the dimen- 
sionless energy, and the radial potential is given by 

if r < R, 
Vm^r) -<(,_, ,. (4c) 

if r> R, 

iTics — s + msgn(eB), (4d) 

where s = R^/{2P) = ^/^o is the magnetic flux 
$ = BR^TT (in units of the magnetic flux quantum 
$0 = h/ |e|) missing inside the circle of radius R. The 
function sgn{-} stands for the sign function. In the nu- 
merical results presented in this paper, we always assume 
that the particle is an electron moving in a magnetic field 
along the positive z-axis, i.e., sgn(eB) = — 1. However, 
our theoretical results are not restricted in such a way. 

Introducing the new variable ^ = r^/2 and transform- 
ing the wave functions in the magnetic region (r > R) as 



/™(r)=Cl™°''^l/'e-«/2;^„(0, (5) 

Eq. 14a|) results in a Kummer differential equation^ 






(1- 



rrics 



<) 






rrics 



-TOcff 



-Xn 



0. 



(6) 



Thus the ansatz for the radial wave function in the mag- 
netic region can finally be written as 



/™(r)=el'"^«-|/'e-«/2 

'1+ I rricif I -rrioff ■ 



xU 



,1+ I mcff |,C 



(7) 



where U is the confluent hypergeometric function^. 
Note that the function U tends to zero as r — > oo. 

It is easy to show that the radial wave function inside 
the circle of radius R (where the magnetic field is zero) 
satisfies the Bessel differential equationil.. Thus the ra- 
dial wave function is given by 



i{t) = JmiV^r), 



(8) 



where Jm{x) is the Bessel function of order ni. 

Matching the radial wave functions inside and outside 
the circle gives a secular equation whose solutions are 
the eigenvalues of the system. The matching conditions 
a,i r ~ R yield 



-y- \ngm{T)\ 

dr -^ ^ ' I r=R/l 



Tr'^^-^^t^n,- (^) 



For a given m this secular equation depends only on the 
dimensionless missing flux s. 



III. SEMICLASSICAL APPROXIMATION: THE 
BERRY- TABOR APPROACH 

We now turn to the semiclassical treatment of the sys- 
tem. Generally, in d dimensions, a system is integrable if 
there are d independent constants of the motion. Usually 
this is the result of the separability of the Hamiltonian: 
In a suitably chosen coordinate system the Hamiltonian 
depends only on separate functions <j)i[qi,pi) of the co- 
ordinates and the conjugate momenta. This means that 
the dynamics can be viewed as a collection of indepen- 
dent one dimensional dynamical systems. The function 
4'{QiiPi) plays the role of the Hamiltonian in each sub- 
system. The one dimensional semiclassical quantization 
procedure can be carried out in each subsystem sepa- 
rately 



1 

2^ 



Ptdqi 



("■ + f)' 



0,1,2,, 



(10) 



where li is the action variable and Ui is the Maslov in- 
dex (for details see, eg, Ref . l40|) . The Maslov index is 
the sum of the Maslov indices of the turning points of 
the classical motion. Smooth or "soft" classical turn- 
ing points (zeros oi pi{qi)) contribute -1-1 to the Maslov 
index, while "hard" classical turning points (inflnite po- 
tential walls) contribute +2. Equation 1)10(1 . the Bohr- 
Sommerfeld quantization condition, is widely used to ap- 
proximate the energy levels of classically integrable sys- 
tems. 

Alternatively, from Eq. (|10|l a semiclassical trace for- 
mula known as the Berry- Tabor formula^^ can be derived 
for the oscillating part of the density of states. For two- 
dimensional systems, this formula can be written as 



d{E) = doiE) 



+ °o COS 

■EE — 

p j=\ IT a 






ij'^p 



3/2 / j(»2,p)^ 



rp2 



dI7 



(11) 



(for the detailed derivation of this expression see Ap- 
pendix 0. Here dQ{E) is the average density of states. 
The p-summation runs over the primitive periodic orbits 
of the system, the j-summation runs over their repeti- 
tions; Sp, Tp and Vp denote the classical action, the time 
and the Maslov-index of orbit p, respectively; n2^p is the 
number of cycles in the motion projected to the action 
variable I2 under one cycle of the orbit; and /i = 5(12, E) 
denotes the action variable /i as a function of the energy 
and l2- 



A. Classical dynamics of the system 

It is easy to show that the classical Hamiltonian in 
polar coordinates (r, ip) , inside and outside the non- 
magnetic region is: 



H = ■^ + V(r) 



(12) 



where Pr and p^ are the canonically conjugate momenta, 
and the radial potential V{r) = !^ Vm{T) is the same as 
in l(^ with the following replacements: 



"^cff 



Pv 
h ' 

s + ^sgn(cB). 
n 



(13a) 
(13b) 



Note that here m and TOoS are continuous classical vari- 
ables. As we shall see below in the semiclassical approxi- 
mation, the canonical momentum p^p is quantized accord- 
ing to the Bohr-Sommerfeld quantization rules H1(J|) . 

Since the Hamiltonian does not depend explicitly on 
if (the system is rotationally invariant), the conjugate 
momentum p^ is a constant of motion. Thus the angular 
action variable becomes 



1 

2^ 



p^dip = p^. 



(14) 



The conjugate momentum inside the anti-dot is in fact 
the angular momentum. Outside the anti-dot there is 
an additional term due to the non-zero vector potential. 
From if = dH/dpip one finds 



The derivation in the frame of classical mechanics is out- 
lined in Appendix ^ Following Ref. |4^ the relation be- 
tween these quantities and the corresponding quantum 
states of the system can be derived from quantum me- 
chanics. It turns out that the same relations hold for the 
cyclotron radius and the guiding center provided p,^ is 
quantized as p,^ = Tim, where m now is an integer. Then, 
rrics is the same as that defined by Eq. Ij4d|l . The same 
results were found by Sim et ali2^. Using Eq. (|17|l we 
shall discuss in detail the correspondence between classi- 
cal orbits and quantum states in Sees. IIII Dl and Hvl 

The turning points given in Eq. (|16bp can be expressed 
in terms of the cyclotron radius and the guiding center: 



I r™* 



/r^* 



j c~ g, if TOoff > 0, 
I g - c, if nicff < 0, 

c+ g. 



(18a) 
(18b) 



We can now classify the classical orbits according to 
the relation between the values of the turning points 
given by (|15|l and the corresponding radius of the cir- 
cular non-magnetic region (in units of /) 



Mr^ip, 



if r<R, 



Pv = 



Mr^ip + eB^:^^, if r > R. 



(15) 



/.. 



in Eq. (jTHJ. To 
= g{I^,E) one 



We now choose /i = Ir and I2 — j-ip 
calculate the radial action variable Ir 
needs to perform the integral of pr — \/ E — V{r) be- 
tween the classical turning points of the radial potential. 
For a given E these turning points can be obtained from 
V{r) = E. Using the same dimensionless variables as in 
Sec. ^ we have one turning point for the potential inside 
the non-magnetic circle: 



V^' 



(16a) 



and for the potential valid outside the circle there are 
two turning points: 



2{e + TOcff) T 2 ^£(2 + 2meff) 



y/e + nicfi + Irrioffl T \/e + Woff - \mcs\, (16b) 



R 



2s. 



(19) 



There are two different cases listed in Tabled For orbits 
of type A the particle outside the non-magnetic region 
moves along a cyclotron orbit and then passes through 
the magnetic field free region as a free particle. In the 
case of orbits of type B the particle does not penetrate 
the non-magnetic region. In this case one can further 
distinguish two additional types of cyclotron orbits de- 
pending on the sign of rrics. The condition r°^* > tr 
listed in Table U and Eq. I|18a(l imply that c— g > R ior 
'rricS > 0, and g — c>R for mcs < 0. From a simple geo- 
metrical consideration it follows that in the first case the 
cyclotron orbits (denoted by Bi) lie outside the circle of 
radius R, while in the latter case the orbits (denoted by 
B2) completely encircle the non-magnetic region. These 
conditions can be rewritten as 



msgn{eB) > V2se > 0, for i?i, (20a) 

-msgn(eS) > \/27e > 0, for B2. (20b) 



where the upper/lower sign of =p distinguishes the first 
and second turning points. Note that r°"* < t™*, and 
the turning points are real if either Wcff > or e > 
— 2moff for rrios < 0. 

For a given energy E and momentum p^p one can cal- 
culate the cyclotron radius g and the guiding center c: 



iV^, 



l\/e + 2m, 



■off- 



(17a) 
(17b) 



For both types Bi and B2, e < m^/2s is valid. In the 
case of orbits of type A we have e > rr? jls. 

We now turn to the calculation of the radial action 
variable Ir- Using the radial potential given by l|lc|l inside 
the non-magnetic circle we find 



ei„(£,r) = - prdr ^ / y^e - V{t) dr 



y 



£T^ — m^ — m arccos 



TWe 



(21a) 



and similarly outside, we have 

0out(e,T) = T Prdr ^ / y/e - V{t) dr 



2 V I 2 



1. , X . / £ + TOoff-T2/2 

- (e + moff) arcsm , , ^=r- 

2 \ y/e [e + 2mofr 



-- |meff I arcsm „ . , , ^ =^ -(aib) 

2 y T^V^(^ + 2"^off) / 

The radial action variables for the orbits of types A and 
B can be expressed in terms of the functions Ojn and 
Oout, and are listed in Table |l| 



TABLE L Classification of the orbits and the corresponding 
radial action variables. See also the text. 



Case 




conditions 




TT jr 


A 


T"o" < TR and 

out ^ 
T2 > T^R, 




eout(e,T2°"*)-eout(£,rR) 
+ ein(e,rR)-ein(e,r^") 


B 




out ^ 
Ti > TR 




eout(e,rr')- Oout (e,rr') 


Note that 

©out 
©out 




-J (£ + mcff-|moff|), (22a) 
-0out(e, rr'). (22b) 



Thus, for orbits of type B, the radial action variable /^ 
can be simplified to 



Ir 



- (e + TOeff 



\mcs\) . 



(23) 



It is clear from (I14II and Tabled that for fixed s, the 
radial action variable Ir is a function of the rescaled en- 
ergy £ and the angular action variable I^ through 771 and 
TOoff . Then, for orbits of type A, the partial derivative in 
the denominator of Eq. Hll|) has a rather simple form 






1 



"^cff 



7r^/2 



se - 



2 moff 



(24) 



However, the amplitude for orbits of type B in Eq. l(TT|l 
cannot be calculated using the second partial derivative 
of /r-, therefore the contribution from these orbits to the 
semiclassical level density is calculated separately in Ap- 
pendix 10 

Knowing the explicit e and I^ dependence of the ra- 
dial action variable 7^, the Bohr- Sommerf eld quantiza- 
tion conditions given by Eq. (|10|1 for orbits of type A can 
be rewritten as 



hm. 



Ir = h\ n+ - 



(25a) 
(25b) 




FIG. 2; Two examples for the basic orbit segments (an arc 
in the magnetic field followed by a straight line inside the 
anti-dot) of orbit type A, and the related angles a, (3 and 7. 
A particle traveling along the segments moves anti-clockwise 
with respect to the center of the anti-dot in the case drawn on 
the left, while in the other example it moves clockwise on the 
straight line inside the anti-dot. Thus in this case the angle 
7 is negative. 



where n = 0, 1, 2, • • • and m is an integer, and the energy- 
dependent radial action variable Ir is given in Table|ll(the 
Maslov indices are z^;p = and f^ = 2, for details see, eg, 
Ref. 40) . Using l|23(l for orbits of type B, the semiclassical 
quantization conditions can be simplified and the energy 
levels are 



E„ 



hiVc 



iTOcffl -TOcff + 1 



(25c) 



where m^g ^ s + msgn^eB), and m and n are integers. 
These levels coincide with the familiar Landau levels in 
a homogeneous magnetic field but the quantum number 
m is replaced by nics. Below in Sec. IlllCl we shall com- 
pare the exact energy levels with those obtained from 
the Bohr- Sommerf eld quantization conditions for orbits 
of types A and B. 



B. Periodic orbits 

To apply the Berry- Tabor formula (|ll|l . one needs to 
describe the possible periodic orbits of the magnetic anti- 
dot system. The periodic orbits of type A can be charac- 
terized by their winding number w (the number of turns 
around the center under one cycle) and the number of 
identical orbit segments rig the orbit can be split up to. 
These segments consist of a circular path outside the anti- 
dot followed by a straight line inside. We introduce the 
angles a,/3 and 7 to characterize these basic orbit seg- 
ments as shown in Fig. |21 These angles always fulfill 



nsiP + l) 
sin(a) 
sin(/?) 



wn, 

R 

p 



(26a) 
(26b) 



(a and (3 are always positive and the sign of 7 follows 
the sign of w). The relations between the indices w,ng 
and the angles characterizing the basic orbit segment are 
summarized in Table Hll for the four possible sub-classes. 
When either w is negative (orbits of type Ai), or the cy- 
clotron radius p is smaller than the radius of the anti-dot 
(orbits of type A2), the angles a, (3 and 7 are fully deter- 
mined by w and ris, since in these cases (3 is definitely 
smaller than 7r/2. On the other hand, when w > Q and 
p > R, one must also specify whether (3 is smaller (orbits 
of type ^3) or larger (orbits of type A4) than 7r/2 to fully 
determine the periodic orbit. 



TABLE II: The different sub-classes of orbits of type A and 
the corresponding relations between the angles defining the 
basic orbit segment and the indices w and Us- Every orbit 
with a negative winding number w falls into sub-class Ai, 
regardless of whether p is smaller or larger than R. The sub- 
class A2 consists of those orbits which have a positive winding 
number and p < R. In case oi w > 0, p > R and /3 < n/2, 
the orbit is of type A3, otherwise it is of type A4. The angle 
a can be obtained from w and Us directly. Then /3 can be 
calculated from a, p and R, and finally 7 from a and /3. 



Ai : w <0 




A2 



w > 0, 
p<R 







P = f3o = arcsin(-^ sing) 



J — a — -K — 13 



l3 = -K-(5o 



^ = a — j3 



The action of periodic orbit p can generally be ex- 
pressed as 



Sp = hkLp + eBAp, 



(27) 



where k = \/2ME/h is the wave number, Lp is the length 
of the orbit and Ap is the area inside the magnetic field. 
In our case 



2ns [pa + i?sin(|7|)] , 



Ap ^ ns< p 



-R' 



a + sgn (2q — tt) 



/? + sgn(2/3-^) 



sin 2a 
2 
sin 2/3 



(28) 



(29) 



Finally, to use (|li|l and H24() , we also need the time period 
r„ and /,„ associated to the orbit, which can be written 



"-'■p 



Tp = 



Lp 2n 

V LUc 

Iv,p/^ = ±%/2secos(7). 



(l + V2^sin(|7|)), (30b) 
(30c) 



where v denotes the velocity and in the latter expression 
the upper sign is for the orbits with ly > and the lower 
sign is for the orbits with w < 0. This expression for I^ 
can be obtained from (|14|l : the angular action variable is 
equal to Pip, and, as already mentioned, inside the anti- 
dot p^ is equivalent to the angular momentum. 



C. Results 

In this section we compare the numerically exact 
energy levels with those calculated from the Bohr- 
Sommerfeld quantization conditions. Similarly, we 
present results for the density of states obtained from 
the Berry- Tabor formula (|11|) . 

The numerically exact energy levels of the magnetic 
anti-dot system are calculated from the secular equation 
for fixed m. Solving Eq. H25|) for e we obtain the en- 
ergy levels in the Bohr-Sommerfeld approximation. The 
results for a given magnetic field are shown in Fig.O The 
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FIG. 3: Exact (crosses) and semiclassical (-1- signs) energy 
levels (in units of TujJc/2) of the circular magnetic anti-dot 
obtained from Eqs. and l|25^ as functions of m for s = 5. 
We take sgn(eB) = -1 as in Ref.l2l 



Therefore the action in our units can be written as 
Sp/n = ns |2£a + 2V2se sin(|7|) 



-fsgn(ei3) 



e a + sgn (2a — tt) 



sin 2a 



2s { (3 + sgn {2(3 - tt) 



sin 2/3 



(30a) 



agreement between the exact and the semiclassically cal- 
culated energy levels is excellent. Our results also agree 
with those presented in Ref. i22i. For large \m\ the energy 
levels tend to the Landau levels, while in the opposite 
case a substantial deviation can be seen. In the latter 
case, the energy levels result from the quantization of or- 
bits of type A. In the work by Sim et ali2^ these states 
were called magnetic edge states. One can see that even 



the low-lying energy levels of these magnetic states can 
be accurately calculated in the Bohr-Sommerfeld approx- 
imation. However, a significant deviation of the eigenval- 
ues of these states from the bulk Landau levels can be 
seen in the figure. The lowest energy level of the mag- 
netic anti-dot system is the state m — Q and n = 0. Note 
that the spectrum can be calculated much more easily 
in the semiclassical approximation than from the exact 
secular equation involving the confluent hypergeometric 
function U. 

Increasing the magnetic field, we experienced slight de- 
viations. These discrepancies may be explained qual- 
itatively in the following way. As the magnetic field 
tends to infinity, the charged particle spends less and 
less time outside the circular region, and in the limiting 
case its motion is described by an elastic reflection from 
the boundary of the magnetic and non-magnetic regions. 
The radial potential becomes a hard wall at r = i?. Thus, 
one of the classical turning points for orbits of type A be- 
comes a hard one and the corresponding contribution to 
the Maslov index tends to 2. Blaschke and Brack'*-'^ ob- 
served a similar situation in circular magnetic billiards. 
Their numerical investigations have conflrmed the argu- 
ment presented above. Here we do not discuss this issue 
further. 

We now present results for the density of states calcu- 
lated from the Berry- Tabor formula Hll|) for two magnetic 
fields given by the missing flux quanta s — 5 and s = 10. 
To evaluate the semiclassical density of states in practice, 
we have regularized the periodic orbit sum in (|11|) with a 
Gaussian smoothing by multiplying the amplitude of the 
orbits with e~^^'' , (where S is infinitesimal), as discussed 
jjj4o^^ This factor suppresses the contribution from the 
long orbits and broadens the delta functions at semiclas- 
sical energies. Substituting Eqs. (|^ and (jSUl) into the 
regularized version of 1)111) . we obtain the semiclassical 
density shown in Fig. 0] plotted together with the numer- 
ically obtained quantum energy levels. The agreement 
between the two results is good for the majority of the 
levels, however, in the case of missing flux quanta s = 10, 
apparent discrepancies can be observed for eaxample at 
energies close to e « 3. We think that a better agreement 
can be obtained for stronger magnetic fields by taking 
into account the magnetic field dependent Maslov index. 
The work along this line is in progress. 



D. Phase diagram, the classical-quantum 
correspondence 

In this section we classify the exact energy levels in 
terms of the classical orbits. In Sec. IIIlAl orbits of types 
A and B have been introduced according to the positions 
of the turning points compared to the magnetic anti-dot. 
Orbits of type A can further be classified into sub-classes 
A1-A4, as it has been shown in Sec. IIII Bl However, in 
that section only periodic orbits have been studied. We 
now go beyond the condition l|2()a|l for periodic orbits. 
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FIG. 4: The quantum mechanical energy levels as a function 
of s, and the semiclassical level density obtained from the 
Berry- Tabor approach as function of e are plotted together. 
The smoothing parameter 5 for the periodic orbit sum was 
S = 0.002, and the summation of the orbits runs from Us = 2 
to Tis^niax = 240 and from lo = to lUmax ~ 120. (Numerical 
experience showed that with these maximum values ns^max 
and «;max the broadened delta functions are quite prominent 
at semiclassical energies ). 



Then the four types of sub-classes can be characterized 
by the angles a, (3 and 7, the cyclotron radius g and the 
guiding center c of the classical orbits. Using simple geo- 
metrical arguments, these angles can be calculated from 
g and c. Using Eq. H17(l (which is valid in the quantum 
case, too) the above classical parameters classifying the 
orbits can be directly related to the quantum states given 
by the quantum number m, n (and so the energy eigen- 
value E„i^n) and the missing fiux quanta s. Thus, the 
conditions given in the first row of Table ^ for the dif- 
ferent types of orbits can be reformulated in terms of the 
particle energy e and the quantum number m. The re- 
sults are summarized in Tablelllll A similar classification 



TABLE III: Conditions in terms of e and m for different sub- 
classes of orbits of types A and B. Here we take sgn(e_B) = 
-1. 



Ai 


A2 


A3 


A4 


Bl 


S2 


£ > m^/2s 


£ < mV2s 


m < 


m > 


m < 


m > s 


£ < 2s 


£ > 2s 


m< 2s 


m > 2s 


£ > 2(m - 


s) 



has been made for electronic states of a circular ring in 
magnetic field'*'^ and for ring-shaped Andreev billiards*^. 
It may also be useful to present the conditions listed 
in Table |n] graphically. In e-m space, the classically al- 
lowed regions corresponding to the different orbits look 
like a 'phase diagram'. In Fig. |S1 such a phase diagram is 
plotted in the space of e and m In s-m space, for a given 
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magnetic field. This phase diagram should be compared 
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FIG. 5; Phase diagram in the space of e and m: the classically 
allowed regions defined by the conditions in Table UTTl for the 
different orbits. Here s — 5 and sgn(e_B) = — f . 

with Fig. 13 the plot of the energy levels from the quan- 
tum mechanical calculation. In this way, the different ex- 
act levels can be classified in terms of the corresponding 
classical orbits. We should stress again that these orbits 
are not necessarily periodic. Examples will be shown in 
Sec.HVI 



IV. COMPARISON OF THE CURRENT 

DISTRIBUTIONS AND THE CLASSICAL 

TRAJECTORIES 

An apparent correspondence between the classical or- 
bits and the quantum states can also be made by calcu- 
lating the current flow patterns in the magnetic anti-dot 
system. The particle (probability) current density^^i^ in 
magnetic field is given by 

J = ^ (*grad** - **gradvl') - ^ A|vI'|^ (31) 

Using the vector potential ^ the current density (in our 
units) for states '^m,n can be written as j — jtp—, where 

J^ir) = ^ \^m,nf [2m - sgn{eB) (r^ - r^) 6 (r - tr)] , 

(32) 
and e^ is the unit vector in the Lp direction. 

Solving the secular equation ^ and then determining 
the normalized eigenstates, the related current densities 
can be calculated from (|32|l . Figures FfilSI 1101 andll llshow 
the current flow patterns for given eigenstates and the 
corresponding classical trajectories of the particle. The 
missing flux quanta is s = 5 and sgn(ei3) = — 1 in all flg- 
ures. In these flgures the current density j(r) at point r is 
represented by an arrow with length proportional to the 
magnitude of the current density and the midpoint of the 



arrow is at the point r. At a given energy and quantum 
number m corresponding to the classical canonical mo- 
mentum p^ — hm the cyclotron radius and the guiding 
center of the classical orbit is calculated from Eq. H17|l . 
Hence the classical trajectory of the particle can be de- 
termined and are shown in the flgures (scaling is in units 
of the magnetic length /). In Figs. 151 and [T^ the radial 
dependence of the current density is plotted for the corre- 
sponding eigenstates. States {m,n) = (—1,0), (0,0) and 
(1, 0) were called magnetic edge states in Ref. |2^ 

In Fig. ini the current inside the non-magnetic region 
flows clockwise (the magnitude of the current density is 
negative in accordance with Fig. O, while outside the 
magnetic dot it flows counterclockwise. The classical tra- 
jectories inside the dot form a 'caustic' and the current is 
enhanced here in the clockwise direction. From Table UTTl 
we flnd that the orbit is of type Ai. The trajectory ap- 
parently also satisfles the conditions given in Table ^ 
(without the requirement for the periodicity of the or- 
bits). 

In Fig. [3 for state to = and n = the current is 
zero inside the magnetic dot. This can easily be seen 
from Eq. (|32|l . while the direction of the current flow is 
counterclockwise outside. The orbit is a limiting case 
between types Ai and A2. Figure |H1 for state m — 1 and 
n — shows a counterclockwise current flow both inside 
and outside the magnetic dot. The current jip is positive 
everywhere as it can be also seen in Fig. |21 The orbit is 
again of type A2 , in accordance with the conditions given 
in Tables in and Uni 




FIG. 6: The current flow (in units of h/{2M)) for state m = 
— 1 and 71 = 0. The classical orbit is of type Ai. Here and 
in Figs. l9ll'^ the missing flux quanta is s = 5, and we choose 
sgn(e_B) = —1. 



One can observe that the current density j^p is not dif- 
ferentiable at r = i?. This is because the theta function 
in Eq. H32|l . Physically, this is a consequence of the step 
function behaviour of the magnetic field. Nevertheless, 
the divergence of the current density vector j is still zero 




FIG. 7: The current flow (in units of fi/{2M)) for state m = 
and n = G. The classical orbit is of type A2. 




FIG. 8: The current flow (in units of h/{2M)) for state m = 1 
and n = 0. The classical orbit is of type A2. 



everywhere. 

Finally, orbits of types A^ and A^ are shown in Figs. 1101 
and [TTl for states {m,n) = (1,6) and (14,2), respec- 
tively. Similarly, for these states the corresponding cur- 
rent densities as functions of the distance from the origin 
are plotted in Fig. E| For both states the current is 
very small inside the magnetic dot. In the case of state 
(rn,n) = (1,6), the trajectories almost cross the origin 
(to is small), while for state (to, n) = (14, 2) only a small 
portion of the trajectory penetrates into the magnetic 
anti-dot regions. The qualitative agreement between the 
current flow patterns and the classical trajectories is, 
again, clearly visible. 

Note that in all of these figures the classical trajec- 
tories are not periodic orbits with energy correspond- 
ing to the given eigenstate. This is not surprising, since 
the quantization should not be applied to periodic orbits 




FIG. 9: The current densities jip as functions of r (in units of 
I) for states shown in Figs. I6I8I The vertical line is at r = R. 




FIG. 10: The current flow (in units of h/{2M)) for state m 
1 and n = 6. The classical orbit is of type A3. 



in real space but to the motion on the two-dimensional 
torus parametrized by the action variables and their 
canonically conjugate angle variables (for details see, eg, 
Ref . l4(J) . In fact, the Berry- Tabor formula (|ll|l suggests 
that an infinite number of periodic orbits with proper 
weights can only result in the correct quantum mechani- 
cal density of states. 

It has been shown by Halperini^ that the probability 
current can be related to the derivative of the energy 
levels with respect to the angular momentum quantum 
number: 



je,„dr 



1 dE„ 



(33) 



/Q h dm 

In Fig. 1131 the integral of the current densities given in 
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FIG. 11: The current flow (in units of h/{2M)) for state m — 
14 and n = 2. The classical orbit is of type A4. 




FIG. 12: The current densities j^ as functions of r (in units 
of I) for states shown in Figs. llUlllI The parameters are the 
same as in Fig. |S| The vertical line is at r = R. 



Eq. H33|) is plotted as functions of m for n = 0, 1, 2. One 
can see that for ?ti <C the total current is zero (these 
are the orbits of type Bi in our classification) and that 
for TO ^ 2s it tends to a constant value (type B2). The 
current is negative for states corresponding to orbits of 
type Ai , while it increases monotonically for states corre- 
sponding to orbits of types A2 — A4. This is in accordance 
with the right hand side of Eq. H33|) using the energy dis- 
persion (the m dependence of the energy levels) plotted 
in Fig. 01 




FIG. 13: The total current Im,n (in units of Ti/(2M)) obtained 
from the integral of j^ given in 1133^ as functions of m for 
n = 0, 1, 2. The parameters are the same as in Fig. |S| 



V. CONCLUSIONS 

In this paper we investigated the energy spectrum of 
the circular magnetic anti-dot systems obtained from 
exact quantum and semiclassical calculations. In the 
proper dimensionless variables the only relevant param- 
eter of the system is the missing flux quantum. The 
system is separable, and in the quantum calculation the 
energy levels are the solutions of the secular equation 
derived from the matching conditions of the wave func- 
tions inside and outside the dot. In the semiclassical 
treatment we presented two different methods. On the 
one hand, the density of states was calculated using the 
Berry- Tabor formula. On the other hand, the energy 
levels obtained from the Bohr-Sommerfeld quantization 
rules. The main difference between the two methods is 
that in the first case one needs to characterize the possi- 
ble periodic orbits in real space, while in the latter, the 
motion of the particle is on a torus in the space of the 
action variables. In our numerical results we compared 
the quantum energy spectrum and that obtained in the 
semiclassical approach. We showed that the energy lev- 
els of the magnetic anti-dot systems obtained from the 
two semiclassical methods were in good agreement with 
the numerically exact quantum results for weak magnetic 
fields. However, by increasing the magnetic field, a slight 
deviation between the exact and the semiclassically ap- 
proximated energy levels can be observed. We argued 
that the reason for this discrepancy may be traced back 
to the fact that the Maslov index should be magnetic 
field dependent. A thorough investigation of this prob- 
lem might be an interesting future work. 

A classification of the energy spectrum for arbitrary 
magnetic fields was presented in terms of the classical 
orbits defined by their cyclotron radius and guiding cen- 
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ter. Such identifications are based on the exphcit re- 
lations between these classical parameters of the orbits 
and the quantum states. The correspondence between 
the quantum states and the classical trajectories can be 
made transparent by drawing a phase diagram with re- 
gions corresponding to six different types of orbits in the 
space of energy and angular momentum quantum num- 
ber. 

Finally, we calculated the current flow patterns for 
eigenstates that correspond to orbits with trajectories 
penetrating into the field-free region. The related classi- 
cal trajectories were also shown for the sake of comparis- 
son. From these results one can see the close correspon- 
dence between the structure of the trajectories and the 
distribution of the current densities obtained from the 
quantum calculations. 

From the energy spectrum of the magnetic anti-dot 
systems one can determine the free energy. The good 
agreement between the semiclassical and quantum treat- 
ment of the system allows us to use semiclassical methods 
in the weak field limit for calculating the energy spec- 
trum. Therefore, the semiclassical approach provides a 
useful starting point for successive studies of thermo- 
dynamic properties, such as magnetization. Moreover, 
the semiclassical approximation can be an effective tool 
for investigating arbitrarily shaped magnetic anti-dots 
(which would be a very difficult task in the quantum 
case) or systems with more complicated magnetic field 
profiles. 
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APPENDIX A: THE BERRY- TABOR FORMULA 

The quantized energies can be recovered if we express 
the Hamiltonian in terms of li 

E{ni,n2, ■■■,nd) = H{Ii,l2, ...Jd) = 
H {h (m + j^i/4) ,ti{n2 + 1^2/4:) , ..., h {ud + Vdl^)) . 

(Al) 

The semiclassical density of states is the density of these 
energies: 

d{E)^ Y^ S{E~E{ni,n2,...,nd)). (A2) 

ni,n2,...,nd=0 



The density of states can be rewritten via the Poisson 
resummation technique 



d{E) 



J d''l6{E-H{h,l2,-,Id))x 

d +OC 

d'^Idt 



i—1 ni — — oo 



E 



X 



„ 27r?i''+i 

ITll .1712 ,...,7nd— — 00 



(A3) 



Here, we used the Fourier expansion of the delta spike 
train. The term rrii = {i = 1, 2, ..., d) can be evaluated 
directly and yields the non-oscillatory average density of 
states. Other terms can be evaluated by the saddle point 
method, when h ^ 0. The saddlepoint conditions select 
the periodic orbits of the system, and the result of the 
integration is 



4-oc 



d{E) = doiE)+J2J2 



(2^) 



(d-l)/2 



2(xp-i)n^d+i)/2 



p J=l 

cos{jSp{E)/h^^jiyp + ^{d-l)) 
^{jTpr-^i-defDp) 



(A4) 



Here p is the index of the primitive periodic orbits, j is 
the number of repetitions, Sp is the classical action along 
the orbit, Tp is the time period of the orbit, and Vp is the 
Maslov index. The quantity Xp is the number of action 
variables of the periodic orbit whose saddle point value is 
zero (Ife = 0), since in this case the Gaussian saddle point 
integral is only one-sided, and its contribution is 1/2 of 
the full Gaussian integral. The matrix Dp is related to 
the second derivative matrix 



det D = det 



dlidlj 
dH{h,...,Ia) 

a/, 



dH(h,...Ja) 
dh 





(A5) 



Equation (|A4() is the generic form of the semiclassical 
density of states in terms of periodic orbits, known as 
the Berry- Tabor formulai^^. 

In two dimensions, very often the Hamiltonian cannot 
be expressed with the action variables explicitly, only the 
implicit function 



l2=g{h,H), 



(A6) 



is available. In this case it is more useful to express the 
quantities in the Berry- Tabor trace formula in terms of 
the derivatives of g. Taking the partial derivative of (jA6p 
with respect to /i yields 



= 



dg{h,H) , dg{h,H)dH{h,l2) 



dh OH 

while the partial derivative of 
gives 



(A7) 



dg{h,H)dH{h,l2 



dH 



dh ' 
6)) with respect to h 



(A8) 



dh 
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The frequencies can be expressed from these equations 
as 



UJl 



UJ2 



dH{h,h) 
dh 

dH{h,h) 
dh 



dh 

dgjIuH) 

dH 



9g{Ii,H) 
dH 



(A9) 
(AlO) 



Periodic orbits are recovered from wi = ^^^ and uj2 = 
^^Y^ . The action /i for a periodic orbit at energy E can 
be obtained by solving equation 



UJl 
UJ2 



Til 

n2 



n2.p 



dh ' 



(All) 



where we introduced rii — jni^p and n2 = jn2,p corre- 
sponding to the primitive orbit. Then the period can be 
expressed simply as 



27m 



2,p- 



dgih.H) 



dH 



(A12) 



The main determinant to be calculated reads 



detD 



d^Hiluh) d^Hiluh) dH{h.l2) 




dl'i dlidh dli 




d^H{Iij2) d^H(Ii.l2) dH(h,l2) 


_ 


dlidl2 dl'i a/2 




dH(h,l2) dH{hJ2) n 
dh dh ^ 






d^H fdHV ^ d^H dHdH 
dll \dh) ' "dhdhdhdh 


d^H (dHy\ , 


5/| \dh) ]■ 


V 



(A13) 



Now, the second derivatives of H can be expressed with 
the second derivatives of g by taking further partial 
derivatives of HA7|I and (|A8I) with respect to h and h- 
Then we can express the second derivatives as 



d'^H 
'dl[ 



d^H 
dhdh 

d^H 



1 



d'^g dg dg d^g f dg 



im 



dHdhdhdH dH^ \dh 



dll [dH) ) ' 
1 / d^g dg d^g dg 



ag_Y KdH^dh dhdH dH 

dH 



I 



d'g 



dll 



m'''^'' 



(A14) 
(A15) 

(A16) 



Using these expressions, the determinant becomes 
1 d^g _ (27rn2,p)3 d^g 



detD 



m 



^dii 



T3 dll 



(A17) 



The density of states in two dimensions is then 



d{E) = do{E) + 

sr^P^ COS {jSp{E)/h 



ij^p 



(d-1)) 



.(A18) 



APPENDIX B: DERIVATION OF THE 

CYCLOTRON RADIUS AND THE GUIDING 

CENTER 

The cyclotron radius can be determined from the en- 
ergy E of the particle. The energy is conserved, and 
obviously E — ^ MloI g^ , thus 



(Bl) 



The guiding center may be calculated as follows. As we 
have seen, the conjugate momentum given by Eq. H15|l is 
a constant of motion, therefore, e.g., for r > R the right 
hand side of ((T5|l should also be a constant at any point 
of the orbit. At first apply this equation for points P and 
Q, which are the points closest to and farthest from the 
origin (the center of the circle of radius R) of an orbit 
lying outside the anti-dot. These are special points of 
the orbits for which the right hand side of Eq. I|15|) has 
a simpler form. Then the distances of points P and Q 
from the origin are rp = c— g and rq = c+ g (we assume 
that point Q is farther from the origin). From a simple 
geometrical argument one finds that the angular velocity 
at points P and Q satisfies the following equations 



(B2) 
(B3) 



rpifp = guJcSgn{eB), 
rqipQ = -guJcSgn-{eB). 



Substituting, for example, rp = c— g and ip'p from IJB2I 
into Eq. 115() . and using (|Bip . we find 



- = \/e + 2mcs- 



(B4) 



The same results can be obtained by using l|B3|l for point 
Q. If the oribit encompasses the anti-dot then the right 
hand side of ljB2p should be multiplied by a factor of —1. 
The case of orbits with trajectories penetrating into the 
anti-dot can be treated similarly. However, the expres- 
sions for the cyclotron radius and the guiding center are 
the same as above for all cases. 



APPENDIX C: CONTRIBUTION OF THE 

CYCLOTRON ORBITS TO THE 
SEMICLASSICAL DENSITY OF STATES 

In the case of the cyclotron orbits, the integral in /^ 
in Eq. IjASp has to be calculated directly rather than 
using the saddlepoint method. As I^ is constant, the 
integrand does not depend on the integration variable 
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and therefore the integral is equal to the measure of the 
interval of the possible I^s. Without loss of generality, 
we take sgn(eB) = — 1 in this section. 



and their Maslov index is /it = 2, therefore the contribu- 
tion to the semiclassical level density from these orbits 
reads 



1. Cyclotron orbits of type Bi 

Cyclotron orbits which do not encompass the anti-dot 
(type Bi) are possible at any value of p and any neg- 
ative angular momentum p^p (see Table llllll . At point 
P of these orbits (points P and Q of a cyclotron orbit 
are defined in Appendix IB|| . from Eqs. (|14|) and l|15|l we 
obtain 



lip = Pip = Mrpifip + eB- 



K' 



(CI) 



This is minimal when rp — R (and the cyclotron orbit 
touches the boundary of the anti-dot), and is maximal 
when the orbit is placed as far as possible from the anti- 
dot. By denoting the radius of the system with £ (in 
units of I) the intergation in ljA3p with respect to I^p 
yields a factor A/^^ = I,p{R) — Iip{C). Using ljCl|l and 
(lB2ll. one finds 



A/,. 



1 



= C^/e - V2^e + -{C^ - 2s) 



(C2) 



The Ir- and i-integrals can be evaluated with the saddle 
point method, just as in the case of a one-dimensional 
system. The determinant of the second derivative matrix 
is 



d- {E) = A~ J2 cos (J^^ - J^) ■ (C8) 

i=i 
The sum has Dirac delta peaks at e = 2n + 1, where 

n = 0, 1, 2, . . . and m < 0. These are the familiar Landau 
levels of an electron for m < 0. 



2. Cycltoron orbits of type B2 

For cyclotron orbits encompassing the anti-dot (type 
B2) the angular momentum satisfies the condition 
Pip/h — m > 2s (see Table IITT)l . At point Q of these 
orbits, using (fT^ and ifTsll . we can write 



lip = Pip ^ MrgipQ + eB- 



K' 



(C9) 



The minimum and the maximum of tq are g and 2q— R, 
respectively. Between these values, I^p, as a function of 
rq, is monotonic, thus the integration in IjASp over I^p 
gives A/^ = I^{g) - I^{2g - R). Using (El and (lB3t . 
we have 



A/,„ 



^2S£. 



(CIO) 



d-'H 



detD = det 



rp ati 



dH \ 



~dTf^ dl,. \ 

.m. / 



— — ) . (C3) Similarly to (|C6|) . the amplitude of the orbits becomes 



From Eqs. (|X8|) and ^^, we find 



Al 



A/^ 1 



+ s- 



/2se 



(Cll) 



dH 

detD = 



dg_ 
OE 



(C4) 
(C5) 



Thus, the total amplitude of these orbits in the periodic 
orbit sum is 

A- . ^'^ ' 



fi Tuo, 
1 



huJr 



Cy/e - V2^ + -{C^ - 2s) 



(C6) 



The action can be calculated from Eq. 1)23(1 . and for 
Pip/h = TO < we have 



S = 2'Klr = fl-KE, 



(C7) 



Using H23|) . the action for Ptp/h = m > 2s is 

S = 2nlr ^h7r{e + 2s- 2m) , (C12) 

Finally, the contribution to the periodic orbit sum of 
these orbits is 

00 
d+{E) ^A+Y^ cos [TTJie + 2s~ 2m) - jn] .(C13) 

The sum has Dirac delta peaks at £ = 2{m — s) -I- 2n -I- 1, 
where m and n are non-negative integers, and m > 2s. 
These are again the familiar Landau levels of an electron 
for TO > 2s. 
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